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Summary. We consider dynamical tunneling between two symmetry-related reg- 
ular islands that are separated in phase space by a chaotic sea. Such tunneling 
processes are dominantly governed by nonlinear resonances, which induce a cou- 
pling mechanism between "regular" quantum states within and "chaotic" states 
outside the islands. By means of a random matrix ansatz for the chaotic part of 
the Hamiltonian, one can show that the corresponding coupling matrix element di- 
rectly determines the level splitting between the symmetric and the antisymmetric 
eigenstates of the pair of islands. We show in detail how this matrix element can 
be expressed in terms of elementary classical quantities that are associated with 
the resonance. The validity of this theory is demonstrated with the kicked Harper 
model. 



1.1 Introduction 

Since the early days of quantum mechanics, tunneling has been recognized 
as one of the hallmarks of the wave character of microscopic physics. The 
possibility of a quantum particle to penetrate an energetic barrier represents 
certainly one of the most spectacular implications of quantum theory and has 
lead to various applications in atomic and molecular physics as well as in 
mesoscopic science. Typical scenarios in which tunneling manifests are the es- 
cape of a quantum particle from a quasi-bound region, the transition between 
two or more symmetry-related, but classically disconnected wells (which we 
shall focus on in the following), as well as scattering or transport through 
potential barriers. The spectrum of scenarios becomes even richer when the 
concept of tunneling is generalized to any kind of classically forbidden tran- 
sitions in phase space, i.e. to transitions that are not necessarily inhibited 
by static potential barriers but by some other constraints of the underlying 
classical dynamics (such as integrals of motion). Such "dynamical tunneling" 
processes arise frequently in molecular systems [1] and were recently realized 
with cold atoms propagating in periodically modulated optical lattices [2,3]. 
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Despite its genuinely quantal nature, tunneling is strongly influenced, if 
not entirely governed, by the structure of the underlying classical phase space 
(see Ref. [4] for a review). This is best illustrated within the textbook example 
of a one-dimensional symmetric double-well potential. In this simple case, 
the eigenvalue problem can be straightforwardly solved with the standard 
Wentzel-Kramers-Brillouin (WKB) ansatz [5]. The eigenstates of this system 
are, below the barrier height, obtained by the symmetric and antisymmetric 
linear combination of the local "quasi- modes" (i.e., of the wave functions 
that are semiclassically constructed on the quantized orbits within each well, 
without taking into account the classically forbidden coupling between the 
wells) , and the splitting of their energies is given by an expression of the form 



Here E is the mean energy of the doublet, V{x) represents the double well 
potential, m is the mass of the particle, u denotes the oscillation frequency 
within each well, and the integral in the exponent is performed over the whole 
classically forbidden domain, i.e. between the inner turning points of the or- 
bits in the two wells. Preparing the initial state as one of the quasi-modes 
(i.e., as the even or odd superposition of the symmetric and the antisymmet- 
ric eigenstate), the system will undergo Rabi oscillations between the wells 
with the frequency AE/h. The "tunneling rate" of this system is therefore 
given by the splitting (1.1) and decreases, keeping all classical parameters 
fixed, exponentially with 1/h, what gives rise to the statement that tunneling 
"vanishes" in the classical limit. 

The above expression for the splitting can also be derived in a geomet- 
ric way which is independent of the particular representation of the phase 
space. For this purpose, it is necessary to realize that the two symmetric 
wells are connected in the complexified classical phase space. This is most 
conveniently expressed in terms of the local action-angle variables (I, 9) of 
the well: describing the quantized torus with the action variable I = I n 
in the left well by (pz,(/ n , 6), <7l(/„, 9)) and its counterpart in the right 
well by (pn(I n ,9),qB(I n} 9)) (p and q are the position and momentum vari- 
ables), one can show that the analytic continuations of the two manifolds 
(PL/R(In,9),qL/n(I n ,9)) coincide when 9 is permitted to assume complex 
values [6]. Generalizing standard semiclassical theory [7] to this complex La- 
grangian manifold allows one to reproduce Eq. (1.1), where the exponent now 
contains the imaginary part of the action integral J pdq along a path that 
connects the two tori in complex phase space [8]. 

This approach can be generalized to multidimensional, even non separable 
systems, as long as their classical dynamics is still integrable [8]. It breaks 
down, however, as soon as a non integrable perturbation is added to the 
system (e.g. if the one-dimensional double-well potential is exposed to a pe- 
riodically time-dependent driving). In that case, invariant tori may, for weak 
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perturbations, still exist due to the Kolmogorov-ArnoPd-Moser (KAM) the- 
orem. It can be shown, however, that their analytic continuation to complex 
(multidimensional) angle domain generally encounters a natural boundary in 
form of weak singularities that arise at a given value of the imaginary part of 
the angle variable [9]. Only in very exceptional cases, the complex manifolds 
that originate from the two symmetry-related tori happen to meet in form of 
an intersection [10]. In such a situation, one can apply a semiclassical method 
introduced by Wilkinson, which leads to an expression of the type (1.1) for 
the splitting with a different ft-dependence in the prefactor [11]. 

Indeed, numerical calculations of model systems in the early nineties 
clearly indicated that tunneling in non integrable systems is qualitatively dif- 
ferent from the above one-dimensional case. If the perturbation introduces an 
appreciable chaotic layer around the separatrix between the two wells (which 
in the Poincare surface of section might still be located far away from the quan- 
tized tori), the tunnel splittings generally become strongly enhanced compared 
to the integrable limit. Moreover, they do no longer follow a smooth exponen- 
tial scaling with 1/ft as expressed by Eq. (1.1), but display huge, quasi-erratic 
fluctuations at variations of ft or any other parameter of the system [12, 13]. 

These phenomena are traced back to the specific role that chaotic states 
play in such systems [14-17]. In contrast to the integrable case, the tunnel 
doublets of the localized quasi-modes are, in a mixed regular-chaotic system, 
no longer isolated in the spectrum, but resonantly interact with states that 
are associated with the chaotic part of phase space. Due to their delocalized 
nature, such chaotic states typically exhibit a significant overlap with the 
boundary regions of both regular wells. They may therefore provide an efficient 
coupling mechanism between the quasi-modes - which becomes particularly 
effective whenever one of the chaotic levels is shifted exactly on resonance with 
the tunnel doublet. This coupling mechanism generally enhances the tunneling 
rate, but may accidentally also lead to a complete suppression thereof, arising 
at very specific values of ft or other parameters [18, 19]. 

The validity of this "chaos-assisted" tunneling picture was essentially con- 
firmed by a simple statistical ansatz in which the quantum dynamics within 
the chaotic part of the phase space was represented by a random matrix from 
the Gaussian orthogonal ensemble (GOE) [14,15,20]. In presence of small 
coupling coefficients between the regular states and the chaotic domain, this 
random matrix ansatz yields a truncated Cauchy distribution for the proba- 
bility density to obtain a level splitting of the size AE. Such a distribution is 
indeed encountered in the exact quantum splittings, which was demonstrated 
for the two-dimensional quartic oscillator [20] as well as, later on, for the 
driven pendulum Hamiltonian that describes the tunneling process of cold 
atoms in periodically modulated optical lattices [21,22]. The random matrix 
ansatz can be straightforwardly generalized to the tunncling-induced decay of 
quasi-bound states in open systems, which is relevant for the ionization pro- 
cess of non dispersive electronic wave packets in resonantly driven hydrogen 
atoms [23]. Chaos-assisted tunneling is, furthermore, not restricted to quan- 



4 Peter Schlagheck 1 , Christopher Eltschka 1 , and Denis Ullmo 2 ' 3 

turn mechanics, but arises also in the electromagnetic context, as was shown 
by experiments on optical cavities [24] and microwave billiards [25,26]. 

A quantitative prediction of the average tunneling rate was not possible in 
the above-mentioned theoretical works. As we shall see later on, this average 
tunneling rate is directly connected to the coupling matrix element between 
the regular and the chaotic states, and the strength of this matrix element 
was unknown and introduced in an ad-hoc way. A natural way to tackle this 
problem would be to base the semiclassical description of the chaos-assisted 
tunneling process on complex trajectories [27], to be obtained, e.g., by solving 
Hamilton's equations of motion along complex time paths. Such trajectories 
would possibly involve a classically forbidden escape out of (and re-entrance 
into) the regular islands as well as classically allowed propagation within the 
chaotic sea. 

For the case of time-dependent propagation processes, such as the evolu- 
tion of a wave packet that was initially confined to a regular region, this ambi- 
tious semiclassical program can be carried out in a comparatively straightfor- 
ward way, which is nevertheless hard to implement in practice [28,29]. Indeed, 
Shudo, Ikeda and coworkers showed in this context that the selection of com- 
plex paths that contribute to the semiclassical propagator requires a careful 
consideration of the Stokes phenomenon, in order to avoid "forbidden" trajec- 
tories that would lead to an exponential increase (instead of decrease) of the 
tunneling amplitude [28,29]. This semiclassical method can be generalized to 
scattering problems in presence of non integrable barriers [30-32] and allows 
one to interpret structures in the tunneling tail of the wave function (such as 
plateaus and "cliffs" ) in terms of chaos in the complex classical domain. 

A crucial step towards the semiclassical treatment of "time-independent" 
tunneling problems, such as the determination of the level splitting between 
nearly degenerate states in classically disconnected wells, was undertaken by 
Creagh and Whelan. They showed that the splitting-weighted density of states 
J2 n AE n 5(E — E n ) (where E n are the mean energies and AE n the splittings 
of the doublets) can be expressed as a Gutzwiller-like trace formula involv- 
ing complex periodic orbits that connect the two wells [33,34]. Such orbits 
also permit to determine the individual splittings AE n , provided the wave 
functions of the associated quasimodes are known [35]. In practice, this ap- 
proach can be successfully applied to fully chaotic wells that are separated 
by an energetic barrier [33,34,36], since in such systems the semiclassical 
tunneling process is typically dominated by a single instanton-type orbit. A 
generalization to chaos-assisted tunneling, where many different orbits would 
be expected to contribute, does not seem straightforward. 

In time-dependent propagation problems as well as the in approaches based 
on the Gutzwiller trace formula [37], the Lagrangian manifolds that need to 
be constructed generally have simple analytical structures. This is the case for 
time-dependent problems because the initial wave function is usually chosen as 
a Gaussian wave packet or a plane wave, which are semiclassically associated 
with simple manifolds. In Gutzwillcr-likc approaches, this due to the inherent 
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structure of the theory, but at the cost of involving long orbits when individual 
eigenlevels need to be described. For the problem of determining individual 
level splittings in chaos-assisted tunneling, however, the natural objects with 
which one would have to start with are the invariant tori of the system, the 
analytical structure of which is extremely involved in case of a mixed regular- 
chaotic dynamics. As was mentioned above, these tori cannot be analytically 
continued very far in the complex phase space as a natural boundary is usually 
encountered. As a consequence, it does not seem that a description in terms 
of complex dynamics could be easily developed for chaos-assisted tunneling. 

The main purpose of our contribution is to show that a relatively com- 
prehensive description of the tunneling between two regular island in a mixed 
phase space can nevertheless be obtained. This description involves a num- 
ber of approximations which might be improved to reach better accuracy. It 
however leads to quantitative predictions for the tunneling rates which are in 
sufficiently good agreement with the exact quantum data to ensure that the 
mechanism that underlies this process is correctly accounted for. We shall, 
in this context, particularly focus on the classically forbidden transition from 
the regular island into the chaotic sea. As already pointed out above, the 
associated coupling matrix element determines the average tunnel splittings 
between the islands, which means that a simple semiclassical access to this 
matrix element would open the possibility to quantitatively estimate tunnel- 
ing rates in systems with mixed dynamics. 

An important first step in this direction was undertaken by Podolskiy and 
Narimanov: By assuming a perfectly clean, harmonic-oscillator like dynamics 
within the regular island and a structureless chaotic sea outside the outermost 
invariant torus of the island, a semiclassical expression of the form 

was derived for the average eigenphase splitting [38]. Here A is the phase 
space area covered by the regular island, and r(a, x) denotes the incomplete 
Gamma function [39]. The prefactor 7 is system specific, but does not depend 
on h, which permits a prediction of the general decay behavior of the splittings 
with 1/h. Good agreement was indeed found in a comparison with the exact 
splittings between near-degenerate optical modes that are associated with a 
pair of symmetric regular islands in a non integrable micro-cavity [38] (see 
also Ref. [40]). 

The theory was furthermore applied to the dynamical tunneling process 
in periodically modulated optical lattices, for which the splittings between 
the left- and the right-moving stable eigenmodes were calculated in Ref. [21]. 
Those splittings seem to be very well described by Eq. (1.2) for low and 
moderate values of 1/h, but display significant deviations from this expression 
deeper in the semiclassical regime [38] . Indeed, the critical value of 1/h beyond 
which this disagreement occurs coincides with the edge of an intermediate 
plateau in the splittings, which extends over the rather large, "macroscopic" 
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range 10 < 1/h < 30. As we shall see in Section 1.3, such prominent plateau 
structures appear quite commonly in chaos-assisted tunneling processes (see 
also Ref. [41]), and a smooth, monotonously decreasing expression of the type 
(1.2) cannot account for them. 

To understand the origin of such plateaus, it is instructive to step back 
to the conceptually simpler case of nearly integrable dynamics, where the 
perturbation from the integrable Hamiltonian is sufficiently small such that 
macroscopically large chaotic layers are not yet developed in the Poincare 
surface of section. In such systems, the main effect of the perturbation consists 
in the manifestation of chain-like substructures in the phase space, which 
arise at nonlinear resonances between the eigenmodes of the unperturbed 
Hamiltonian, or, in periodically driven systems, between the external driving 
and the unperturbed oscillation within the well. In a similar way as for the 
quantum pendulum Hamiltonian, such resonances induce additional tunneling 
paths in the phase space, which lead to couplings between states that are 
located in the same well [42,43]. 

The relevance of this effect for the near-integrable tunneling process be- 
tween two symmetry-related wells was first pointed out by Bond et al. [44] 
who argued that such resonances may lead to a strong enhancement of the 
tunneling rate, due to couplings between lowly and highly excited states 
within the well which are permitted by near-degeneracies in the spectrum. 
In Refs. [45,46], a quantitative semiclassical theory of near-integrable tun- 
neling was formulated on the basis of this principal mechanism. This theory 
allows one to reproduce the exact quantum splittings on the basis of purely 
classical quantities that can be extracted from the phase space, and takes 
into account high-order effects such as the coupling via a sequence of different 
resonance chains [45,46]. Recent studies by Keshavamurthy on classically for- 
bidden coupling processes in model Hamiltonians that mimic the dynamics of 
simple molecules confirm that the "resonance-assisted" tunneling scenario pre- 
vails not only in one-dimensional systems that are subject to a periodic driving 
(such as the "kicked Harper" model which was studied in Ref. [45,46]), but also 
in autonomous systems with two and even three degrees of freedom [47,48]. 

Our main focus in this review is that such nonlinear resonances play an 
equally important role also in the mixed-regular chaotic case. Indeed, they 
have recently been shown to be primarily responsible for the coupling between 
the regular island and the chaotic sea in the semiclassical regime [49]. In 
combination with the above random matrix ansatz for the chaotic states, a 
simple analytical expression for the average tunneling rate is obtained in this 
way, which provides a straightforward interpretation of the plateau structures 
in the splittings in terms of multi-step coupling processes induced by the 
resonances. 

To explain this issue in more detail, we start, in Section 1.2, with a descrip- 
tion of the effect of nonlinear resonances on the quantum dynamics within a 
regular region in phase space. We then discuss how the presence of such reso- 
nances leads to a modification of the tunnel coupling in the near-integrable as 
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well as in the mixed regular-chaotic regime, and present a simple semiclassical 
scheme that allows one to reproduce the associated tunneling rates. Applica- 
tions to tunneling processes in the kicked Harper model are studied in Section 
1.3. 



1.2 Theory of resonance-assisted tunneling 

We restrict our study to systems with one degree of freedom that evolve under 
a periodically time-dependent Hamiltonian H(p,q,t) = H(p,q,t + r). We 
suppose that, for a suitable choice of internal parameters, the classical phase 
space of H is mixed regular-chaotic and exhibits two symmetry-related regular 
islands that are embedded into the chaotic sea. This phase space structure is 
most conveniently visualized by a stroboscopic Poincare section, where p and 
q are plotted at the times t — nr(n e Z). Such a Poincare section typically 
reveals the presence of chain-like substructures within the regular islands, 
which arise due to nonlinear resonances between the external driving and 
the internal oscillation around the island's center. We shall assume now that 
the two islands exhibit a prominent r:s resonance — i.e., where s internal 
oscillation periods match r driving periods, and r sub-islands are visible in 
the stroboscopic section. 

The classical motion in the vicinity of the r:s resonance is approximately 
integrated by secular perturbation theory [50] (see also Ref. [46]). For this 
purpose, we formally introduce a time-independent Hamiltonian Ho(p, q) that 
approximately reproduces the regular motion in the islands and preserves 
the discrete symmetry of H . The phase space generated by this integrablc 
Hamiltonian consequently exhibits two symmetric wells that are separated 
by an energetic barrier and "embed" the two islands of H. In terms of the 
action-angle variables (I, 9) describing the dynamics within each of the wells, 
the total Hamiltonian can be written as 



H(i,e,t) = H (i) + v(i,e,t) 



(1.3) 



where V would represent a weak perturbation in the center of the island. 

The nonlinear r:s resonance occurs at the action variable I r:s that satisfies 
the condition 

2tt .,, _ dff 



f f ■ Q S ~ 



with Q r:s = 



dl 



(1.4) 



i=i r 



We now perform a canonical transformation to the frame that corotates with 
this resonance. This is done by leaving / invariant and modifying 9 according 
to 



d 



(1.5) 



This time-dependent shift is accompanied by the transformation H ^ H = 
H — Q r:s I in order to ensure that the new corotating angle variable $ is 
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conjugate to /. The motion of I and d is therefore described by the new 
Hamiltonian 

H{I,$,t)=H {I)+V{I,$,t) (1.6) 

with 

H (I) = H (I) - Q r ..J , (1.7) 

V{I,&,t) = V(I,#+n r:e t,t). (1.8) 

The expansion of Ho in powers of / — I r:s yields 

H (I) * H^ + (/ ~ /r " )2 + O ((/ - /, :s ) 3 ) (1.9) 

with a constant H ^ and a quadratic term that is characterized by the effective 
"mass" parameter m r . s . Hence, dHo/dl is comparatively small for / ~ I r:s , 
which implies that the corotating angle d varies slowly in time near the res- 
onance. This justifies the application of adiabatic perturbation theory [50], 
which effectively amounts, in first order, to replacing V(I,$,t) by its time 
average over r periods of the driving (using the fact that V is periodic in t 
with the period rr) [51]. By making a Fourier series expansion for V(I, 9, t) in 
both 9 and t, one can show that the resulting time-independent perturbation 
term is (27r/r)-periodic in i9. We therefore obtain, after this transformation, 
the time- independent Hamiltonian Ho(I) + V av (I, i?) where V av can be written 
as the Fourier series 

1 f rT °° 
V av (I,d) = — V(J,0,t) = yV fc (I)cos(fcr0 + fc )- (1.10) 
rT Jo k=0 

This effective Hamiltonian is further simplified by neglecting the action de- 
pendence of the Fourier coefficients of V av — i.e., we use Vk = Vk(I = I r -.s) in 
Eq. (1.10) — and by employing the quadratic approximation (1.9) of H (I) 
around I = I r:s . Leaving out constant terms, we finally obtain the effective 
integrable Hamiltonian 

(I - T ) 2 °° 

H eS (1, 0) = K r:s > +J2 V k cos ( fc ^ + • (I-") 
Zm r . s fe=i 

The structure of this Hamiltonian exhibits all ingredients that are neces- 
sary to understand the enhancement of tunneling. This can be most explicitly 
seen by applying quantum perturbation theory to the semiclassical quantiza- 
tion of H c ff, treating the "kinetic" term cx ( I — I r:s ) 2 as unperturbed part 
(I = —ihd/d-d) and the "potential" term Y^k=i ^ cos(kri) + fa) as perturba- 
tion: Within the unperturbed eigenbasis ($\n) ~ exp(\nd) (i.e., the eigenbasis 
of H in action- angle variables), couplings are introduced between the states 
|n) and |n + kr) according to 
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{n + kr\H eS \n) = ^V k e^ . (1.12) 

As a consequence, the "true" eigenstates \4> n ) of H e s contain admixtures from 
unperturbed modes \n') that satisfy the selection rule \n'— n\ = kr with integer 
k. They are given by the perturbative expression 



hM = |n)+£ ' \n + kr) + 

k^o n n+kr 



E 



k k^O ^ n E'n+kr E n E n -\-kr-\-k f r 



\n + kr + k'r) + ... (1.13) 



where E n denote the unperturbed eigenenergies of H c g, i.e., the eigenvalues 
of H Q (I) - Q r . a I. 

Within the quadratic approximation of Hq(I) around I r:s , we obtain from 
Eq. (1.11) 

E « = (L ") 
where the quantized actions are given by 

I n = h(n + 1/2) (1.15) 

(taking into account the generic Maslov index fi = 2 for regular islands). This 
results in the energy differences 

E n - E n , - —L_ (J„ - /„,)(/„ + /„, - 27 r:s ) . (1.16) 

/From this expression, we see that the admixture between \n) and \n') becomes 
particularly strong if the r:s resonance is symmetrically located between the 
two tori that are associated with the actions /„ and /„/ — i.e., if /„ + /„/ ~ 
2I r . s . The presence of a significant nonlinear resonance within a region of 
regular motion provides therefore an efficient mechanism to couple the local 
"ground state" — i.e, the state that is semiclassically localized in the center 
of that region (with action variable Iq < I r - S ) — to a highly excited state 
(with action variable Ik r > Ir-.s)- 

It is instructive to realize that the Fourier coefficients Vk of the perturba- 
tion operator decrease rather rapidly with increasing k. Indeed, one can derive 
under quite general circumstances the asymptotic scaling law 

V k - (fcr)To exp[-fcr/2 r:s i im (/ r:s )] (1.17) 

for large k, which is based on the presence of singularities of the complexi- 
fied tori of the integrable approximation Ho(I) [46]. Here ti m (T) denotes the 
imaginary time that elapses from the (real) torus with action I to the nearest 
singularity in complex phase space, 7 corresponds to the degree of the sin- 
gularity, and V contains information about the corresponding residue near 
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the singularity as well as the strength of the perturbation. The expression 
(1.17) is of little practical relevance as far as the concrete determination of 
the coefficients V k is concerned. It permits, however, to estimate the relative 
importance of different perturbative pathways connecting the states |n) and 
|n + kr) in Eq. (1.13). Comparing e.g. the amplitude A2 associated with a 
single step from \n) to \n + 2r) via V2 and the amplitude Ai associated with 
two steps from \n) to |n + 2r) via V\, we obtain from Eqs. (1.16) and (1.17) 
the ratio 

07^2-7^2 

A 2 /A 1 ~ T e i(^-20O (L18) 
m r:s Vo 

under the assumption that the resonance is symmetrically located in be- 
tween the corresponding two tori. Since V can be assumed to be finite in 
mixed regular-chaotic systems, we infer that the second-order process via the 
stronger coefficient V\ will more dominantly contribute to the coupling be- 
tween \n) and \n + 2r) in the semiclassical limit ft — > 0. 

A similar result is obtained from a comparison of the one-step process 
via Vk with the fc-step process via Vi, where we again find that the latter 
more dominantly contributes to the coupling between \n) and \n + kr) in the 
limit ft — > 0. We therefore conclude that in mixed regular-chaotic systems the 
semiclassical tunneling process can be adequately described by an effective 
pendulum-like Hamiltonian in which the Fourier components Vk with k > 1 
are completely neglected: 

g cff (J, fl) = (/ ~ /r:s)2 + 2V r ..s cosn? (1.19) 

with 2V r:s = V\ [49] (we assume 0i = without loss of generality) . This simple 
form of the effective Hamiltonian allows us to determine the parameters I r:s , 
m r:s and V r:s from the Poincare map of the classical dynamics, without explic- 
itly using the transformation to the action-angle variables of Hq. To this end, 
we numerically calculate the monodromy matrix M r:s = d(pf, qf)/d(pi, m) of 
a stable periodic point of the resonance (which involves r iterations of the 
stroboscopic map) as well as the phase space areas S+ s and S~ s that are en- 
closed by the outer and inner separatrices of the resonance, respectively (see 
also Fig. 1.1). Using the fact that the trace of M r:s as well as the phase space 
areas S^ s remain invariant under the canonical transformation to (1, 1?), we 
infer 

I r:s = ±-(S+ S + S- S ), (1.20) 
\ -'»',: X ,- - ^{S+ a - S-J , (1.21) 
^-arccos(trM r:s /2) (1.22) 



m r . s 



from the integration of the dynamics generated by H e g. 
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Fig. 1.1. Classical phase space of the kicked Harper Hamiltonian (at r = 3), showing 
a regular island with an embedded 9:8 resonance. The thick solid line represents the 
"outer" and the thin dashed line the "inner" separatrix of the resonance. 

Quantum mechanically, the tunneling process between the nth excited 
quantized torus and its counterpart in the symmetry-related island manifests 
itself in a small level splitting between the associated symmetric and antisym- 
metric eigenstates. In our case of a periodically driven system with one degree 
of freedom, these eigenstates arise from a diagonalization of the unitary time 
evolution operator U over one period r of the driving, and the splitting is 
defined by the difference 

A<p n = \<p+ - <p~\ (1.23) 

between the corresponding eigenphases tp^ °f the symmetric and antisym- 
metric state. In the integrable limit, these eigenphase splittings are trivially 
related to the energy splittings AEn^ of the unperturbed Hamiltonian H via 

A^=rAE^/h. (1.24) 

The latter can be semiclassically calculated by the analytic continuation of 
the tori to the complex domain [8] , and are given by 

AE (0) = Mn eM _ an/n) (L25) 
7T 

where fl n is the oscillation frequency of the nth quantized torus and a n de- 
notes the imaginary part of the action integral along the complex path that 
joins the two symmetry-related tori. 

In the near-integrable case [45,46], the presence of a prominent r:s res- 
onance provides an efficient coupling between the ground state and highly 
excited states within the regular region. Within perturbation theory, we ob- 
tain for the modified ground state splitting 
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^o=^ 0) +Ei4 r ; s) i 2 M°, ) (i-26) 
k 

where A^ r '^ = (kr\ip a ) denotes the admixture of the (fcr)th excited component 
to the perturbed ground state according to Eq. (1.13). The rapid decrease of 
the amplitudes with k is compensated by an exponential increase of the 

unperturbed splittings Atp^J , arising from the fact that the tunnel action a n 
in Eq. (1.25) generally decreases with increasing n. The maximal contribution 
to the modified ground state splitting is generally provided by the state \kr) 
for which Ikr+Ia — 2Ir-.s — i-e., which in phase space is most closely located to 
the torus that lies symmetrically on the opposite side of the resonance chain. 
This contribution is particularly enhanced by a small energy denominator (see 
Eq. (1.16)) and typically dominates the sum in Eq. (1.26). 

In the mixed regular-chaotic case, invariant tori exist only up to a maxi- 
mum action variable I c corresponding to the outermost boundary of the reg- 
ular island in phase space. Beyond this outermost invariant torus, multiple 
overlapping resonances provide various couplings and pathways such that un- 
perturbed states in this regime can be assumed to be strongly connected 
to each other. Under such circumstances, the classically forbidden coupling 
between the two symmetric islands does not require any "direct" tunneling 
process of the type (1.25); it can be achieved by a resonance-induced transition 
from the ground state to a state within the chaotic domain. 

The structure of the effective Hamiltonian that describes this coupling 
process is depicted in Fig. 1.2. We assume here that the perturbation induced 
by the nonlinear r:s resonance is adequately described by the simplified pen- 
dulum Hamiltonian (1.19). Separating the Hilbert space into an "even" and 
"odd" subspace with respect to the discrete symmetry of H and eliminating 
intermediate states within the regular island leads to an effective Hamiltonian 
matrix of the form 



/ E V cS ••• \ 
: : 



(1.27) 



V H N1 H NN I 

for each symmetry class. The effective coupling matrix element between the 

fc-i 

Eq — Ei r 



ground state and the chaos block (H^j) is given by 



V cS = V r:s ft V ™ (1.28) 



where E n are the unperturbed energies (1-14) of -ff e ff- Here \kr) represents 
the lowest unperturbed state that is connected by the r:s resonance to the 
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Eq Vr:s 
Vr : s Ef Vr : s 
Vr:s '•• 




V r:3 
Vr:s -B(fc-l)r ' ■ 

Vr:s 
Vr:s E ) 



Fig. 1.2. Sketch of the effective Hamiltonian matrix that describes tunneling be- 
tween the symmetric quasi-modes in the two separate regular islands. The regular 
parts (upper left and lower right band) includes only components that are coupled 
to the island's ground state by the r:s resonance. The chaotic part (central square) 
consists of a full sub-block with equally strong couplings between all basis states 
with actions beyond the outermost invariant torus of the islands. 



ground state and located outside the outermost invariant torus of the island 

(i.e., I(k-l)r < Ic < hr)- 

In the simplest possible approximation, which follows the lines of Refs. [15, 
20] , we neglect the effect of partial barriers in the chaotic part of the phase 
space [14] and assume that the chaos block (H^ ) is adequately modeled by 
a random hermitian matrix from the Gaussian orthogonal ensemble (GOE). 
After a pre-diagonalization of (H^), yielding the eigenstatcs </>^ and eigenen- 
ergies we can perturbatively express the shifts of the symmetric and an- 
tisymmetric ground state energies by 

Performing the random matrix average for the eigenvectors, we obtain 

\{kr\^)\ 2 ^l/N (1.30) 
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for all j = 1 . . . N, which simply expresses the fact that none of the basis 
states is distinguished within the chaotic block (H^j). 

As was shown in Ref. [20], the random matrix average over the eigenval- 



ues £f gives rise to a Cauchy distribution for the shifts of the ground state 



3 



energies, and consequently also for the splittings 

AE = \E+-E ~ 1 (1.31) 

between the symmetric and the antisymmetric ground state energy. For the 
latter, we specifically obtain the probability distribution 



P{AEq) = 1 AE » (1.32) 

1 °> 7T (AE ) 2 + (AE ) 2 

with 

AEi= |f (1.33) 

where A c denotes the mean level spacing in the chaos at energy E . This 
distribution is, strictly speaking, valid only for AE <C V c s and exhibits a 
cutoff at AE ~ 2V c g, which ensures that the statistical expectation value 
(AEo) = J °° xP(x)dx docs not diverge. 

Since tunneling rates and their parametric variations are typically studied 
on a logarithmic scale (i.e., log (AEo) rather than AEo is plotted vs. 1/H, see 
Figs 1.3-1.5 below ), the relevant quantity to be calculated from Eq. (1.32) 
and compared to quantum data is not the mean value (AEo), but rather 
the average of the logarithm of AE . We therefore define our "average" level 
splitting (AEo) g as the geometric mean of AE , i.e. 

(AEo) g = exp [(ln(AEo))] (1.34) 

and obtain as result the scale defined in Eq. (1.33), 



(AE ) S = AE . (1.35) 

This expression further simplifies for our specific case of periodically driven 
systems, where the time evolution operator U is modeled by the dynamics 
under the effective Hamiltonian (1.27) over one period r. In this case, the 
chaotic cigenphases ipj- = S^t/H are uniformly distributed in the interval 
< ipf < 2ir. We therefore obtain 

for the mean level spacing near E . This yields 

<ZW 3 ^(^o) 9 =(^) 2 (1-37) 
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for the geometric mean of the ground state's eigenphase splitting. Note that 
this final result does not depend on how many of the chaotic states do actually 
participate in the sub-block (H^ ); as long as this number is sufficiently large 
to justify the validity of the Cauchy distribution (1.32) (see Ref. [20]), the 
geometric mean of the eigenphase splitting is essentially given by the square 
of the coupling V c s from the ground state to the chaos. 

The distribution (1.32) also permits the calculation of the logarithmic 
variance of the eigenphase splitting: we obtain 

([ln(^ )-(ln(Z\^))] 2 ) = ^. (1-38) 

This universal result predicts that the actual splittings may be enhanced or 
reduced compared to (Aip ) g by factors of the order of exp(7r/2) ~ 4.8, inde- 
pendently of the values of ft and external parameters. Indeed, we shall show 
in the following section that short-range fluctuations of the splittings, arising 
at small variations of ft, are well characterized by the standard deviation that 
is associated with Eq. (1.38). 



1.3 Application to the kicked Harper model 

To demonstrate the validity of our approach, we apply it to the "kicked 
Harper" model [52] 

oo 

H(p,q,t) ~ cos p + ^ rS(t — tit) cosq . (1.39) 



n— — oo 



This model Hamiltonian is characterized by the parameter r > which cor- 
responds to the period of the driving as well as to the strength of the pertur- 
bation from integrability. The classical dynamics of this system is described 
by the map (p, q) ^> (p 1 , q') = T(p, q) with 

p'=p + Tsinq (1-40) 
q' = q — Ts'mp' (1-41) 

that generates the stroboscopic Poincare section at times immediately before 
the kick. The phase space of the kicked Harper is 2ir periodic in position and 
momentum, and exhibits, for not too large perturbation strengths r, a region 
of bounded regular motion centered around (p, q) = (0, 0) (see the upper panel 
of Fig. 1.3). 

The associated time evolution operator of the quantum kicked Harper is 
given by 

U = cxp ( — ^cospj exp ( — ^cos<?j (1-42) 
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where p and q denote the position and momentum operator, respectively. The 
quantum eigenvalue problem drastically simplifies for h = 2tt/N with integer 
N > 0, since the two phase-space translation operators T\ — exp(27rip / 'h) and 
T 2 = cxp(—2niq/h) mutually commute with U and with each other in that 
case [52] . This allows us to make a simultaneous Bloch ansatz in both position 
and momentum — i.e., to choose eigenstates with the properties 

V>(<7 + 2tt) = V(<?) exp(i£,) (1.43) 
+ 2tt) = $(p) exp(i€p) (1.44) 

where ip denotes the Fourier transform of ip. Since the subspace of wave func- 
tions satisfying (1.43,1.44) at fixed Bloch phases £ g and £ p has finite dimension 
N, finite matrices need to be diagonalized to obtain the eigenstates of U. 

Quantum tunneling can take place between the central regular region 
around (0, 0) and its periodically shifted counterparts. The spectral mani- 
festation of this classically forbidden coupling process is a finite bandwidth of 
the eigenphases ip n = (pn q '^ of U that are associated with the nth excited 
quantized torus within this region. We shall not discuss this bandwidth in the 
following (the calculation of which would require diagonalizations for many 
different values of £ q and £ p ), but consider a simpler, related quantity, namely 
the difference 

^ n = |^°)-^°)| (1.45) 

between the eigenphases of the periodic (£ g = 0) and the anti-periodic (£ 9 = it) 
state in position, at fixed Bloch phase £ p = in momentum. In this way, we 
effectively map the tunneling problem to a double well configuration, with the 
two symmetric wells given e.g. by the regions around (0, 0) and (2w, 0). 

Fig. 1.3 shows the eigenphase splittings of the kicked Harper in the near- 
integrable regime at r = 1. The splittings were calculated for the nth excited 
states at N = 6(2n + 1), i.e. for all possible states that are, in phase space, 
localized on the same classical torus with action variable tt/6. These states 
were identified by comparing the overlap matrix elements of the eigenstates of 
U with the nth excited eigenstate (as counted from the center of the region) 
of the time-independent Hamiltonian 

T 

Ho(p, q) = cosp + cosq — — sinpsing 

T 2 T 3 

— — (cospsin 2 q + cosgsin 2 p) — — sin(2p) sin(2g) , (1.46) 
12 48 

which represents a very good integrable approximation of H at small r [46] . 
Multiple precision arithmetics, based on the GMP library [53], was used to 
compute eigenphase splittings below the ordinary machine precision limit. 

The semiclassical calculation of the eigenphase splittings is based on three 
prominent nonlinear resonances that are located between the torus with action 
variable 7r/6 and the separatrix: the 16:2 resonance [54], the 10:1 resonance, 
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Fig. 1.3. Classical phase space (upper panel) and quantum eigenphase splittings 
(lower panel) of the kicked Harper at r = 1. The latter are calculated for the nth 
excited states at N = 2n/h = 6(2n + 1) — i.e., for all states that are semiclassically 
localized on the torus with action variable tt/6. The decay of the exact quantum 
splittings (circles) is quite well reproduced by the semiclassical prediction (thick 
solid line) that is based on the 16:2, the 10:1, and the 14:1 resonance (highlighted in 
the upper panel, together with the above torus) . The dotted and dashed lines show, 
respectively, the semiclassical splittings that are obtained by taking into account 
only the 10:1 resonance, and the "unperturbed" splittings Aip^ that would result 
from the integrable approximation (1.46). 
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and the 14:1 resonance. Below N ~ 100, we find that the tunneling process is 
entirely induced by the 10:1 resonance (i.e., the most dominant one, according 
to the criterion that r and s be minimal). In this regime, the eigenphase split- 
tings are reproduced by an expression of the form (1.26) (with the resonance 
parameters I r:s , m r:Sl V r:s extracted from the classical phase space), where the 
unperturbed splittings Aipn^ are derived, via Eq. (1.25), from the imaginary 
action integrals a n along complex orbits between the two symmetry-related 
regular regions. 

The other two resonances come into play above N ~ 100, where they "as- 
sist" at the transition across the 10:1 resonance. In that regime, a recursive 
application of the resonance-assisted coupling scheme, taking into account all 
possible perturbative pathways that involve those resonances (see Ref. [46]), is 
applied to calculate the semiclassical tunnel splittings. We see that the result 
systematically overestimates the exact quantum splittings above TV ~ 200, 
and does not properly describe their fluctuations. We believe that this mis- 
match might be due to incorrect energy denominators and coupling matrix 
elements that result from the simplified form (1.19) of the effective pendu- 
lum Hamiltonian. The average exponential decay of the quantum splittings, 
however, is well reproduced by the semiclassical theory. A comparison with 
the unperturbed splittings Aipffl calculated from the intcgrable approxima- 
tion (1.46) (dashed line in Fig. 1.3) clearly demonstrates the validity of the 
resonance-assisted tunneling mechanism. 

At r = 2, the phase space of the kicked Harper becomes mixed regular- 
chaotic, and the regular region around (0, 0) turns into an island that is em- 
bedded into the chaotic sea. Fig. 1.4 shows the phase space together with the 
corresponding eigenphase splittings. The latter were calculated here for the 
ground state of the island, with the dimension N ranging from 4 to 300 in 
integer steps. We clearly see that the fluctuations of the splittings are much 
more pronounced than in the near-integrable regime. 

The semiclassical calculation of the eigenphase splittings is based on a 
prominent 8:2 resonance within the island (visible in the upper panel of 
Fig. 1.4) and evaluated according to Eq. (1.37). The sharp steps of (Aip } g 
arise from the artificial separation between perfect regularity inside and per- 
fect chaos outside the island: when ft drops below the value at which the 
(fcr)th excited state of the island is exactly localized on the outermost invari- 
ant torus, (k + 1) instead of k perturbative steps are required to connect the 
ground state to the chaotic domain, and the corresponding coupling matrix 
element V c s (1-28) acquires an additional factor. In reality, the transition to 
the chaos is "blurred" by the presence of high-order nonlinear resonances and 
"Cantori" [55,56] in the vicinity of the island. The latter provide efficient bar- 
riers to the quantum flow [57-59] and therefore change the structure of the 
block (H^ ) in the matrix (1.27). Hence, except for the case of a very "clean", 
structureless chaotic sea, the steps are not expected to appear in this sharp 
form in the actual, quantum splittings. 
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Fig. 1.4. Classical phase space (upper panel) and quantum eigenphase splittings 
(lower panel) of the kicked Harper at r = 2. The splittings are calculated for the 
ground state within the regular island, as a function of N = 2n/h. The thick solid 
line in the lower panel represents the semiclassical prediction of the eigenphase 
splittings, which is based on the 8:2 resonance displayed in the upper panel. The 
size of the logarithmic standard deviation according to Eq. (1.38) is indicated by the 
two dashed lines that accompany the semiclassical splittings; these dashed lines are 
defined by (A<po)g x exp(±7r/2). The long-dashed curve represents the prediction of 
the eigenphase splittings according to the theory of Ref. [38] (see Eq. (1.47)). 
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It is therefore remarkable that the quantum splittings exhibit plateaus at 
approximately the same levels that are predicted by the resonance-assisted 
mechanism. These plateaus, however, seem to be "shifted" to the left-hand 
side with respect to the semiclassical splittings — i.e., the latter apparently 
overestimate the position of the steps. We believe that this is due to a rich 
substructure of partial barriers within the chaos, which effectively increases 
the phase space area of the region in which quantum transport is inhibited. 
Such a substructure is not visible in the Poincare surface of section, but its ex- 
istence becomes indeed apparent when individual trajectories are propagated 
in the vicinity of the island. A more refined approach, taking into account 
such partial barriers in the chaos, would probably be required to obtain a 
better agreement. 

The central fixed point (0,0) becomes unstable at r = 2, when the 2:1 
resonance emerges in the center of the island. This 2:1 resonance develops 
into a symmetric pair of regular islands located along the p = —q axis, which 
dominate the phase space for r > 2.5 together with their counterparts from 
the bifurcation of the fixed point at (7r,7r). The stable periodic points of the 
2:1 resonance are fixed points of the twicely executed kicked Harper map 
T 2 (i.e., where Eq. (1.40,1.41) is applied twice). Therefore, quantum states 
associated with those islands appear, at fixed Bloch phases £ g , £ p , as doublets 
in the eigenphase spectrum of the corresponding time evolution operator U 2 . 
We shall show now that the splitting between the levels of such a doublet is 
again described by the resonance- and chaos-assisted tunneling scenario in the 
semiclassical regime. 

The doublets associated with the ground state of the 2:1 resonance islands 
are calculated by diagonalizing U 2 and by identifying those eigenstates that 
most sharply localized around the centers of the two islands (we restrict our- 
selves to the pair that is located on the q = —p axis). Fig. 1.5 shows the 
corresponding eigenphase splittings as a function of N = h/2n, computed 
for periodic boundary conditions £ 9 = £ p = at r = 2.8 (left column) and 
r = 3 (right column). The semiclassical calculation of the eigenphase split- 
tings is performed in the same way as for r = 2, with the little difference 
that 2r rather than r is used as period in Eq. (1.37). For the coupling to the 
chaos, we identify a relatively large 5:4 sub-resonance within the 2:1 resonance 
islands at r = 2.8, and a smaller 9:8 sub-resonance at r = 3. 

It is instructive to notice that a rather small variation At = 0.2 of the 
perturbation parameter can lead to qualitatively different features of the tun- 
neling rates (compare the lower panels of Fig. 1.5), and that these features are 
indeed reproduced by the resonance-assisted tunneling scheme. At t — 2.8, 
the 5:4 sub-resonance is located sufficiently closely to the center of the island 
that it induces near-resonant internal transitions within the island , i.e. cou- 
plings from the ground state to the (5Z)th excited states that are strongly 
enhanced in the expression (1.28) due to nearly vanishing energy denomina- 
tors Eq — E51 . As a consequence, pronounced peaks in the tunneling rate are 
obtained in the vicinity of such internal near-degeneracies [60] . Such peaks are 




Fig. 1.5. Chaos- assisted tunneling in the kicked Harper at r = 2.8 (left column) 
and t = 3 (right column). Calculated are the eigenphase splittings between the sym- 
metric and the antisymmetric quasimodes that are localized on the two bifurcated 
islands (i.e., the two islands facing each other on the q = —p axis), for periodic 
boundary conditions £ 9 = f p = 0. The semiclassical calculation (thick lines in the 
lower panels) is based on a 5:4 resonance at r = 2.8 and on a 9:8 resonance at 
r = 3. The huge peaks in the lower left panel are induced by internal near-resonant 
transitions within the island, which are permitted since the 5:4 resonance lies rather 
closely to the center of the island. As in Fig. 1.4, the two dashed lines accompanying 
the semiclassical splittings indicate the size of the logarithmic standard deviation 
according to Eq. (1.38), and the long-dashed curve represents prediction according 
to the theory of Ref. [38] . 
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indeed displayed by the exact quantum splittings as well, though not always 
at exactly the same position and with the same height as predicted by semi- 
classics. A completely different scenario is encountered at r = 3, where the 
coupling to the chaos is mediated by a tiny 9:8 sub-resonance that is closely 
located to the outermost invariant torus of the island. This high-order res- 
onance induces rather larges plateaus in the chaos-assisted tunneling rates, 
which are clearly manifest in the quantum splittings, and which indicate that 
the coupling to the chaos is essentially governed by the same matrix element 
over a wide range of ft. 

In addition to the resonance-assisted eigenphase splittings (1.37), we also 
plot in Figs. 1.4 and 1.5 the prediction that is based on the semiclassical 
expression (1.2) proposed by Podolskiy and Narimanov [38]. This expression 
involves the unknown rate 7 which, however, is independent of ft and char- 
acterizes the tunneling process also in the deep anticlassical regime ft ~ 1. 
Hence, we conclude that, for dimensionality reasons, 7 has to be of the order 
of the intrinsic scale t -1 of our system, up to dimensionless prefactors of the 
order of unity. We therefore set 77- = 1 and obtain 

v -r(i/ + i,o)-^3 e {lA<) 

as prediction for the average eigenphase splittings, with v = A/(ttH) where A 
is the phase space area covered by the island. In all cases that were studied 
in this work, we find good agreement of Eq. (1.47) with the exact quantum 
splittings for comparatively low 1/ft, and significant deviations deeper in the 
semiclassical regime. This indicates that the simple harmonic-oscillator ap- 
proximation for the dynamics within the island (which is needed in order 
to predict the tunneling tail of the eigenfunction) is correct for large ft, but 
becomes invalid as soon as nonlinear resonances come into play. 

While the exact quantum splittings may, depending on the size of ft, con- 
siderably deviate from both semiclassical predictions (1-37) and (1-47) of the 
average, the short-range fluctuations of the splittings seem to be well charac- 
terized by the universal expression (1.38) for the logarithmic variance: As is 
shown by the dashed lines accompanying (Z\<po) 9 in Figs. 1.4, the amplitudes 
of those fluctuations are more or less contained within the range that is de- 
fined by the standard deviation associated with Eq. (1.38). This confirms the 
general validity of the chaos-assisted tunneling scenario. 



1.4 Conclusion 

In summary, we have presented a straightforward semiclassical scheme to re- 
produce tunneling rates between symmetry-related regular islands in mixed 
systems. Our approach is based on the presence of a prominent nonlinear res- 
onance which induces a coupling mechanism between regular states within 
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and chaotic states outside the islands. The associated coupling matrix ele- 
ment can be directly extracted from classical quantities that are associated 
with the resonance. Assuming the presence of a structureless chaotic sea, a 
random matrix ansatz can be made for the chaotic part of the Hamiltonian, 
which results in a simple expression for the average tunneling rate in terms of 
the above matrix element. 

Application to the kicked Harper model shows good overall agreement and 
confirms that plateau structures in the quantum splittings originate indeed 
from the influence of nonlinear resonances. A significant overestimation of the 
quantum splittings is generally found for "weakly" chaotic systems where a 
large part of the phase space is covered by regular islands. This is tentatively 
attributed to the presence of a rich substructure of partial barriers (such as 
Cantori or island chains) in the chaotic sea. Such partial barriers generically 
manifest in the immediate vicinity of regular islands and inhibit, for not too 
small h, the quantum transport in a similar way as invariant tori [57,58]. 
Semiclassical studies in the annular billiard [16, 17] do indeed indicate that 
quantum states localized in this hierarchical region [59] play an important role 
in the dynamical tunneling process [61]. 

In the "anticlassical" regime of large h, our semiclassical theory system- 
atically underestimates the exact quantum splittings. This might be due to 
the approximations that are involved in the present implementation of the 
resonance-assisted tunneling scheme (we neglect, e.g., the action dependence 
of the Fourier coefficients of the effective potential in Eq. (1.10), which could 
play an important role in this regime). It is also possible that the coupling 
to the chaos is induced by a different mechanism at large h, which effectively 
amounts to extracting the associated matrix element from the overlap of the 
tunneling tail of the local semiclassical wave function with the chaotic phase 
space. Indeed, we find that the simple semiclassical expression (1.2) intro- 
duced by Podolskiy and Narimanov [38], which is essentially based on that 
scheme, reproduces the tunneling rates quite well in this regime. 

The validity of the resonance-assisted tunneling mechanism was confirmed 
not only for the kicked Harper, but also for the kicked rotor [49] as well as 
for the driven pendulum Hamiltonian that describes dynamical tunneling of 
cold atoms in periodically modulated optical lattices [21]. The theory can be 
furthermore generalized to describe chaos-assisted decay processes of quasi- 
bound states in open systems, such as the ionization of non dispersive wave 
packets in microwave-driven hydrogen [62]. Ongoing studies on dynamical 
tunneling in autonomous model Hamiltonians with two and three degrees 
of freedom [47, 48] clearly reveal that nonlinear resonances play an equally 
important role in more complicated systems as well. The mechanism presented 
here might therefore provide a feasible scheme to predict, understand, and 
possibly also control dynamical tunneling in a variety of physical systems. 
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